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SECTION  I 
INTRODUCTION 


The  efficient  design  of  missiles  requires  accurate  predictions  of 
aerodynamic  characteristics.  .The  early  missile  designers ^were  forced  to 
■test  models  of  theif  design  in  wind  tunnels  to  obtain  these  predictions. 

But  as  the  flow  about  bodies  became  better  described  by  theory,  methods, 
of  calculating  the  aerodynamic  characteristics  were  devised.  However, 
these  methods  are  valid  only  for  a  rather  narrow  ran^e  of  flow  conditions. 
As  the  theories  were  advanced  and  improved',  the  predictions  became  more 
accurate.,  and  the  limitations  with  regard  to  such^parameters  as  angle  of 
attack  and; Mach  number  bacame  less  restrictive.  This  is  not  to  mean  that 
the  Mach  number  and  angle?of-attack  ranges  'Were  large. 

Recently,. .changing  requirements  for  the  purpose  and  general  structure 
of  missiles  require  predicting  the  aerodynamic  characteristics  for  smaller 
1 ength- to-di ameter  ratios  over  larger  ranges  of  angle  of  attack  and  Mach 
'number.  This  study  provides,  advances  to  existing  computational  techniques 
and  furnishes  greater  flexibility  with  regard  to  these  parameters.  The 
results  are  limited  .to  basic  nose-cylir.der  combinations  featuring  pointed 
noses  and  blunt  bases  because  no  evaluations  of  blunt  nose  -  cylinder 
combinations  were  included. 

To  accomplish  the  objective  of  the'  study,  a  computerized  document 
search  provided  about  1J1  ,000  report  abstracts,  and  about  1  ,400  were  select¬ 
ed  for  more  careful  e; /mi nation.  Of  these,  some  400  were  obtained  for 
further  study.  A  review  of  these  reports  showed  that  very  little  experi¬ 
mental  data  was  available  in  the  region  of  interest;  therefore,  the  effort 
for  the  study  was  shifted  from  an  empirical  study  to  a  review  and  more 
thorough  examination  of  theories.  An  analysis  of  the  more  successful 
methods  indicated  that  a  digital  computer  program  .might  be  produced  to 
easily  apply  these  methods  to  simple  model  combinations.  By  avoiding  the 
limiting  assumptions  that  were  commonly  imposed,  it  appeared  that  the 
accuracy  could  be  extended  to  the  region  of  interest. 

This  study,  then,  represents  the  first  stage  of  an  effort  to  produce 
a  method  that  Will  accurately  predict  aerodynamic  characteristics  for 
complex  missile  configurations. 
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SECTION  II 


REVIEW  OF  THEORIES 


1.  MUNK'-S  POTENTIAL  THEORY 

Most  of  the  successful  theories;  related  to  this  study  have  begun  With 
the,  classical  work  of  Mu'.ik  presented  in  Reference  1  and  expanded  by 
Kelly  in  Reference  2  arid  3.  Munk  used  potential  flow  theory  to  show  that 
the  force  distribution  or/ an  airship  iis<:given  by: 


,dF 


(k2  -  ’JCjj  sin  2a.  dx 


0) 


In  this  expression,  S  is  the  surface  area,  a  is  the  anqle  of  attack, 
x  is  the  distance  along  the  longitudinal  axis,  and  (k2  -  ;jc  ;  is  the 
apparent  mass  factor  for  the  difference  in  longitudinal  and  transverse 
flow.  Equation  (1)  integrates  to  zero  for  an  airship,  but  if  the  body 
has  $  flat  base  of  .area  Sg,  the  result  is  a  finite  force: 

Lift  =  jdF  =  lijjv  (k2  -  k^)  Sg  sin  2a  .(2) 


Since  Munk's  research  was  concerned  with  very  small  angles  of  yaw, 
this  potential  force  is  referenced  to  the  wind  vector  as  lift.  At  higher 
angles  of -attack,  this  estimate  would,  yield  increasing  error  because  the 
forces  due  to  viscous  effects  had  not  been  included. 

2.  ALLEN'S  WORK 

Allen,  in  Reference  4,  assumes  that  the  viscous  forces  are  independent 
of  the  potential  flovy,  and  the  lift  and  moment  coefficients  are  written 
as  the  sum  of  two,  terms;  namely,  the  potential  term  and  a  viscous  term. 


S  A 

CL  =  (k2  -  kj  ^  sin  2a  cos  |  +  n  Cd  ^  sjn2a  cosa 

,  _rV-SB^-xm> 


m 


A 


R 


A. 


(k2  -  kx)  sin  2a  cos  | 


+  n  cdc  5^  <xm  '  x>  s1"2n 


(3a) 

(3b) 
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The  first  term  is  Munk's  potential  theory  as  modified  by  Ward  in 
Reference  5  to  show  that  the  potential  cross  force  is  directed  midway 
between  (the  normals  to  the  logitudinal  axis,  and  the  wind  vector.  The  AR 
represents  the  reference  area- for  the  coefficients,  V  is  the  volume  of  K 
the  model,  and  l  is  the  total  length.  In  the  second (viscous)  term,  C, 
refers  to  the  cylinder  steady-state  two-dimensional  drag  ac 

coefficient  for  two-dimensional  flow.  The  value  of  n  corrects  this  drag 
coefficient  for  finite  fnodel  length  according  to  values  tabulated  by 
Goldstein  in  Reference  6.  The  x  is  the  moment  reference  point,  x  is 
the  centroid  of  the  planform  m  area,  A  ,  and  sin2a  comes  from  cthe 
resolution  of  dynamic  pressure  from  the  p  crossflow  velocity  U  sin  a  to 
the  'free-stream  velocity... 

3.  KELLY'S  INTERPRETATION  ■ 

Kelly,  in  References  2  and  3,  deals  with  the  normal  force  coefficient 
which  eliminates  the  factor  cos  a  in  the  second  term  of  Equation  (3a). 

This  method  also  eliminates  the  need  to  assume  that  the  axial  drag  compon¬ 
ent  of  the  lift  force  is  negligible. 

Allen,;  in  Reference  4,  presents  the  viscous  contribution  to  the  cross 
force  from  a  cylindrical  element  of  length  dx  as: 


d.F  =  2r  C.  (p  V  2/2)  dx  (4) 

c 

where  r  is  the  body  radius  at  point  x  (r  =  r(x)),>  and  V  is  the  cross-flow 
velocity.  Kelly  resolved  the  free-stream  dynamic  c  pressure,  q, 
into  cross-flow  and  axial  components  arid  corrected  the  cross-flow  drag 
coefficient  for  finite  length  with  n»  Defining  the  normal  .force  and  moment 
coefficient  as: 


CN  =  N/qA 

(5a) 

Cm  =  M/qAd 

(5b) 

(where  d  is  the  body  diameter),  the  coefficient  increments  due  to  viscous 
cross-flow  are: 

=  <*> 


3 


(6b) 


i 


-  2n  sin2a 

V 


x)>  dx 


These  terms  are  entirely  due  to  viscous  effects,  since  0,  will  not  exist 
for  a  true  potential  flow.  Qc 

Allen  was  aware  that  the  viscous  term  drag  coefficient;  C .  should  be 
related  to  the  transient  effect  as  presented  by  Schwabe  in  ac 
Reference  7.  Schwabe  found  that  the  drag  coefficient  for  a  circular 
cylinder  in  cross-flow  starting  from  rest  varies  according  to  the  'parameter 
VU/r,  where,  t  is  elapsed  time  (Figure  1).  The  steady-state  value  of  C. 
for  Schwabe's  experiment  was  approximately  one.  Kelly  considered" a  Qc 
lamina  of  air,  dx  thick,  mov.ing  along  the  body  axis  at  speed  U  cos  a 
and  across  the  body  at  speed  U  sin  a.  He  postulated  that  the  cross-flow 
senses  an  impulsive  flow  about  a  circular  cylinder  of  radius  r  (Figure  2); 
hence,  the  Schwabe  time  parameter  is  redefined  as 


h  +•  U  sin  a  x  _  x  tan  a  /7\ 

r  '  ~  <r  Cl  cos  a  ~  f  K  1 

To  aid  integration,  KelYy  approximated  Schwabe's  curve  (Figure  1)  with  a 
portion  of  an  arc  tangent  function,  y  =  tan^x,  and  then  approximated  the 
arc  tangent  function  with  a  truncated  Taylor  series.  This  permitted  rapid 
evaluation  of  the  function  without  reference,  to  a  table  of  values. 

Allen  had  used  $  laminar  value  of  1.2  for  C.,  ,  based  on  the  cross- 
flow  Reynolds  number.  In  Reference  8,  Hill  ac  notes  that  Kelly  used 
a  turbulent  value  of  0.35  when  the  free-stream  Reynolds  number  indicated 
the  boundary  layer  to  be  turbulent.  The  observation  riiade  was  that  if  the 
axial  boundary  layer  is  turbulent,  the  cross-flow  layer  will  behave  in  a 
turbulent  manner,  even  though  the  cross-flow  Reynolds  number  is  well 
below  that  associated  with  transitions  (References  3,  8). 

To  evaluate  the  factor  (k?  -  k-,),  Kelly  replaced  it  with  (AB/2)  cos2e, 
where  A  a  1  +  k-j ,  B  =  1  +  k«,  and  tan  0  =  dr/dx  (Figure  3).  This 
approach  has  been  compared  to  the  work  presented  by  Upoon  and  Klikoff  in 
Reference  9. 
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Figure  1.  Schwabe's  Transient  Drag  Coefficient 
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Figure  2. 


Vc  =  ,U  sin  a 
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Model  at  Angle  of  Attack  Schematic' 
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Figure  3.  Plot  of  k-i  and  kg  (Munk) 
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4.  SUPERSONIC  POTENTIAL  FLOW 

Except  for  vehicles  with  slender  noses,  Hunk's  theory  Is  not  accurate 
for  supersonic  flow;  hence,  Kelly  used>  the  second  order  theory  presented 
by  Van  Dyke  In  Reference  10  for  this  flow  region.  This  changes  the 
potential  terns  of  the  coefficient  expressions  to 


(8a) 


(8b) 


5.  BOUNDARY  LAYER  CORRECTIONS 

Kelly,  in  his  calculations,  also  included  model  size  increases  due  to 
thickening  of  the  boundary  layer  as  the  flow  progressed  along  the  body 
length.  These  relatively  small  corrections  (about  5%  for  a  14  caliber 
missile  at  a  =  8  degrees,  and  M  =  2.87,  and  about  4%  for  a  9  caliber  model 
at  a  =  30  degrees  and  M  =  0.26)  are  most  applicable  to  slender  bodies  at 
small  angles  of  attack.  As  the  angle  of  attack  becomes  larger,  the 
distance  a  particle  flows  along  the  body  becomes  smaller  until  at  a  =  ;90 
degrees  there  is  no  flow  along  the  body,  >but  only  cross-flow.  These 
corrections  are  generally  small  and  not  significant  at  relatively  large 
angles  of  attack  for  the  models  used  in  this  study;  therefore,  they  will 
not  be  considered  in  this  report. 

6.  PREVIOUS  RESULTS 

Kelly  was  able  to  achieve  excellent  results  up  to  about  30  degrees 
angle  of  attack  for  subsonic  flow,  and  up  to  10  degrees  for  supersonic 
flow.  With  minor  computational  changes,  he  also  has  produced  reasonably 
good  predictions  up  to  40  or  45  degrees  for  both  cases. 
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SEGTION  III 


DISCUSSION 


The  success  achieved  by  Kelly  shows  that  the  basic  approach  is  sound-* 
and  further  extension  is  possible.  A  simple  modular  digital  computer 
program  allows  for  more  precise  evaluation  of  the  various  parameters.  More, 
specfficaVly,  such  a  program  permits  the  rapid  calculation  of  the  aero-? 
dynamic  characteristics  even  though  the  parameters  may  vary  over  the  -large 
ranges. 

1.  TRANSIENT  EFFECT 

The  transient  evaluated  by  SchWabe's  experiments  limited  Kelly's  method 
to  a  time  parameter  value  of  nine  (x/r  tan  a  =  9),  since  Schwabe's  expert 
ments  did  not  define  the  transient  effect  beyond  this  value.  This  limit 
is  expressed  as  angle  of  attack  versus  model  length  in  Figure  4.  Schwabe's 
limiting  factor  was  the  size  of  his  equipment.  He  showed  a  maximum  C. 
of  2.07  at  the  maximum  time  parameter;  The  data  points  show  a  leveling 
off  of  the  drag  coefficient,  but  no  indication  of  how  it  returns  to  a 
steady-state  value. 

The  first  computations  of  the  computer  program  indicated  a  slight 
over-prediction  of  normal  force  coefficient  near  the  parameter  value  of 
nine  for  subsonic  cases,  the  supersonic  cases  showed  an  extensive  over¬ 
prediction  at  nearly  all  time  parameter  values  in  the  transient  region. 

The  drag  buildup  for  a  cylinder  normal  to  impulsive  flow  is  attributed 
to  the  formation  of  vorticies.  As  movement  begins,  an  upstream  displace¬ 
ment  of  the  fluid  occurs  over  a  broad  front.  The  flow  field  quickly  forms, 
and  begins  to  narrow  as  vorticies  are  developed.  As  velocity  increases, 
the  vortices  begin  to  be  shed,  and  the  wake  narrows  to  near  steady-state 
conditions.  The  narrowing  of  the  wake  results  in  a  cross-flow  drag  reduction. 

Schwabe's  experiments  were  conducted  in  a  water  tank  of  a  length  that 
prevented  his  experiment  from  reaching  steady-state  conditions.  It  is  a 
reasonable  assumption  that  the  upstream  effect  was  influenced  by  the  con¬ 
fines  of  the  tank  length  and  resulted  in  an  erroneous  drag  increase  and 
delayed  the  onset  of  the  reduction  towards  steady-state  flow.  Based  on  a 
knowledge  of  these  limitations,  a  reduction  of  the  transient  maximum  to 
1.97  at  a  time  parameter  of  seven  is  used.  A  linear  decline  to  tlid  steady- 
state  drag  coefficient  at  a  time  Darameter  of  ten  is  used  to  provide  con¬ 
tinuity  of  the  function  (Figure  5).  Since  upstream  effects  provide  a 
major  part  of  the  transient  phenomena,  the  supersonic  case  does  not  utilize 
any  transient  effect. 
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Figure  4.  Limit  of  Schwabe's  Results 
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Figure  5.  Approximation  of  Transient  Effect 


S.chwabe's  experiments  were  conducted  at  a  Reynolds  number  of  580, 
for  which  the  steady-state  drag  coefficient  (Reference  11)  is  about  1.15. 
Schwabe's-  curve  is  therefore  normalized  by  this  value  arid  is  used  as  a 
transient  factor  to  apply  to  the  actual  two-dimensional  drag  coefficient 
for  the  model  and  flow  conditions  being  considered. 

2.  CROSS-FLOW  DRAG-  COEFFICIENT 

The  cross-flow  drag  coefficient  is  a  function  of  the  cross-flow  Mach 
number  and:  Reynolds  number.  To  develop  this  function,  two-dimensional 
data  from- several  sources  are  combined.  In  Reference  IT,  Hoerner  shows 
an  assortment  of  data  at  various  Mach  and  Reynolds  numbers.  His  plots  of 
C .  versus.  iMach  number  agree  wi th  most  other  reports ,  but  they  are  not 
a  continuous  where  the  magnitude  of  the  Reynolds  number  causes  the 
transition1  from  laminar  to  turbulent  flow  to  occur  below  a  Mach  number  of 
approximately  0.55.  Transition  causes  a  flat-bottomed  dip  in  the  drag 
coefficient  (Reference  12)  for  Reynolds  numbers  between  approximately 
300,000  and  700,000.  Above  Mach  0.55,  compressibility  effects  appear 
to  dominate  the  flow  and  suppress  any  Reynolds  number  effects. 

To  Compute  this  dip  in  the  cross -flow  drag  coefficient  curve,  a 
Reynold's  number  is  computed.  The  characteristic  length  used  in  this 
computation  is-  that  which  the  freestream  senses,  and  it  is  expressed  as 
l  =  2r  sin  a  +  x  cos  a.  Freestream  velocity  is  used  in  the  computation, 
based  on  evidence  that  the  cross-flow  boundary  layer  appears  to  assume 
the  character  of  the  axial  boundary  layer. 

The  drag  curves  obtained  from  Reference  11  are  approximated  by 
straight  line  segments  for  the  computer  function.  Another  straight  line 
segment  provides  transition  from  the  laminar  curve  to  the  turbulent  curve 
as  the  Reynolds  number  increases  from  400,000  to  500’, 000  (Figure  6). 

3.  END  EFFECTS 

The’  two-dimensional  drag  coefficient  must  be  corrected  for  finite  , 
cylinder  length.  Goldstein's  data  tabulated  in  Reference  6  gives  an 
overall,  or  average,  ratio  of  CQ  to  Cd  (Figure  7). 

inf 

Since  a  portion  of  the  computer  program  developed  for  this  study  uses 
a  strip  technique  (numerical  integration  by  trapezoid  rule),  the  local 
correction  referenced  to  the  nearest  end,  rathdr  than  the  average  value, 
must  be  used.  The  drag  coefficient  distribution  is  not  clearly  indicated 
in  any  of  the  referehtes  examined  except  in  Gowen  and  Perkin's  report 
(Reference  12).  Efforts  to  curve-fit  these  data  (Figure  8),  led  to  the 
choice  of  the  following  algebraic  function  to  express 

n,  =  l-exp[-(Mc  +  lJyTL]  (9) 
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erkins'  Experimental  Results  (from  Reference  12) 


EL  is  the  distance  from  the  nearest  end  of  the  model,  and  the  quantity 
(M  I  +  1)  arises  from  the  fact  that  the  end  effects  diminish  rapidly  with 
increasing  Mach  number.  This  is  confirmed  by  Penland  in  Reference'  13, 
where  end  effects  have  been  determined  to  be  virtually  non-existarit  for 
Mach  6.86. 

Since  these  end  effect  corrections  apply  to  right  circular  cylinders, 
further  corrections  are  required  for  a  vehicle  with  a  pointed  nose.  As 
the  nose  becomes  more  slender,  the  end  effects  increase  and  the  local  drag 
coefficient  returns  to  that  for  a  right  cylinder.  Thus,  the  drag  reduction 
.Hue  to  end  effects  is  further  reduced  by  the  oblique  flow  around  a  pointed 
nose.  If  the  nose  half-angle  is  near  zero  degrees,  the  flow  is  essentially 
the  same  as,' a  cylinder  of  equal  length!  If  the  nose  half-angle  is  near 
90  degrees,  the  end  effect  is  so  drastically  reduced  that  almost  no  drag 
is  contributed  by  the  nose  portion  of  the  model.  These  observations 
indicate  that  the  end  effect  change  due  to  nose  half-angle  follows  the  form 
of  the  cosine  of  the  nose  half-angle.  Therefore,  when  the  strip  being 
computed  lies  on  the  nose  portion  of  the  model,  the  rr  computed  is  reduced 
by  multiplication  with  the  cosine  of  the  nose  half-angle. 

4.  NOSE  RADIUS 

Computation  of  the  transient  effect,  the  nose  half-angle,  and  the  two- 
dimensional  drag  coefficient  requires  that  the  radius  at , any  axial  point 
be  available.  Since  the  commonly  used  nose  shapes  include  cones,  tangent 
ogives,  secant  ogives,  power  series,  and  parabolics,  a  variable  function 
is  used  to  express  radius  variation  with  axial  distance.  Based  on  a  nose- 
type  code,,  the  applicable  function  is  selected  for  each  computation 
(Figure  9). 

A  tangent  ogive  function  is  used  for  both  ogives,  since  including  a 
function  for  a  secant  ogive  would  be  more  complex  than  necessary  at  this 
stage.  Power  series  and  parabolic  nose  shapes  are  approximated  by  power 
series  expression,  and  cones  are  defined  by  a  tangent  function.  None  of 
the  functions  employed  account  for  any  nose  bluntilng;  hence,  this  analysis 
is  valid  only  for  pointed  noses.  However,  minor  blunting  would  introduce 
only  minor  error  in  the  results. 

5.  HALF-ANGLE 

The  computation  of  the  apparent  mr  factor  AB/2  cos20w  requires  an 
average  value  of  0„  (=tan~1  dr/dx)  over  the  entire  model.  vAlso,  the 
computation  of  n-  Requires  "he  use  of  the  nose  half-angle.  (For  power 
series  and  parabilic  noses,  the  local  half-angle  would  increase  accuracy 
but  this  additional  c Duplication  is  not  included  at  this  stage.)  To 
compute  these  factors,  the  average  nose  half-angle  is  computed  by  ten 
increments  of  Ar/Ax,  each  equal  to  one-tenth  of  the  nose  length.  In  the 
computations,  this  value  is  proportioned,  or  weighted,  according  to  nose 
length  divided  by  the  total  model  length  in  older  to  obtain  an  average 
value  of  6  for  the  model. 
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Figure  9.  Nose  Radius  Computation  Methods 


'6.  DCMDA  AND  DANDA 


The  supersonic  potential  flow  term  requires  the  determination^ of 
values  for  the  derivatives  (dC  /da)  =  0  and  (dCN/da)  =  0  for  the 
particular  model  being  considered.  a  If  these  nvalu£s  are  available 
from  experimental  data  or  other  sources,  they  may  be  read  into  the  computer 
program  on  the  input  data  cards.  If  they  are  not  read  in  as  data,  two 
functions  will  be-  used  to  approximate  them. 


The  approximation- of  Grimmiriger,  Williams,  and  Young  (Reference  14)- 
for  ogives  and  cones  and  afterbodies,  (Figure  10)  is  used  as  a  starting, 
point  for  these  functions.  (For  purposes  of  these  functions,  power  series 
and  parabolic  noses  have  been  grouped  together  with  ogives.)  These  curves 
are  further  approximated  by  straight  line  segments,  with  corrections  for 
Mach  number.  Through  logic  switching,  the  result  produced  is  an  approximate 
function  subroutine.  The  application  is  limited  to  cone  half-length 
(inscribed  cone  for  noses  not  true  cones)  from  five  to  fifteen  degrees. 
Comparisons  with  other  data  (Reference  2)  indicate  reasonable  agreement. 
Since  potential  force  is  a  minor  part  of  the  total  normal  force  in  the 
region  of  interest,  the  approximation  is  adequate. 


The  dC  /da  is  computed  by  multiplying  the  dCj,/da  by  the  moment  arm, 
measured  from  the  model  center  of  pressure  to  the^model  base.  Accuracy  of 
the  results  has  suffered  somewhat  since, no  particular  effort  has  been  made 
to  accurately  predict  the  moment  of  the  center  of  pressure  for  varying 
model  dimensions  or  flight  conditions.  The  inclusion  of  these  factors  would 
significantly  complicate  the  method  of  this  study.  When  additional  experi¬ 
mental  data  are  available,  further  study  of  these  refinements  would  be  in 
order. 
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Figure  10.  Initial  Normal  Force  Slope  (From  Reference  14) 
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wMm  iv 

PROGRAM  DESCRIPTION 


The  computer  program  resulting  from  this  study  is  written  for  tha< 

IBM 360  Model  65  as  installed  in  tho  Louisiana  State  University  Computer 
Research  Center.  The  WATFIV  compiler,  which  gives  complete  syntax  and' 
logic  analysis  with  explicit  error  descriptions,  was  used  to,  insure  an 
accurate  and  rapid  transition  from  one  phase  of  the  program  development 
to  the  next.  Some  of  the  special  capabilities  of  this  Toad-and-go  compiler 
were  utilized;  consequently,  changes-  in  the  program  might  be  required  if 
another  type  of  computer  was  used". 

1.  PROGRAM  INFORMATION 

The  computer  program  information  is  provided  by  a  series  of  comment 
cards  which  define  the  variables,  codes  *  and  terms  used  in  the  program 
(Table  I). 

2.  MAIN  PROGRAM 

The  main  program  controls  the  overall  operation  of  the  computer  program- 
(Table  II).  This  program  provides  a  means  of  data  input,  terminates 
execution,  writes  output  headings  (which  include  input  data),  and  exercises 
logic  to  differentiate  between  subsonic  and  supersonic  computation  techniques. 

The  input  data  for  each  case  are  read  from  a  single  card.  They  include 
a  nose-type  code,  the  nose  length  and  model  length  (both  in  calibers),  the 
freestream  Mach  number,  the  Reynolds  number  per  caliber  divided  by  the  Mach 
number,  and  the  constants  (dC^/da)  =  0  and  (dC  /da)  =  0  if  experimental 
or  other  source  data  are  to  be  use3.  The  program  termination  is  based  on 
a  nose-type  code  of  zero,  which  is  supplied  by  a  blank  card  at  the  end  of 
the  data  stack.  The  next  statements  call  the  proper  subroutine  for  the  case 
being  considered,  and  then  returns  to  the  input  step  for  either  the  next 
case  or  termination. 

3.  SUBSONIC  SUBROUTINE 

The  subroutine  that  directs  the  computation  for  the  subsonic  case  is 
shown  in  Table  III.  It  computes  constants,  initializes  incremental  variables, 
increments  the  angle  of  attack,  computes  the  first  (potential)  term,  calls 
another  subroutine  for  the  second  (viscous)  term,  and  prints  out  the  com¬ 
puted  values.  All  dimensional  constants  are  computed  in  calibers.  Some 
constants  are  computed  through  the  use  of  additional  subprograms,  which  will 
be  described. 
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TABLE  I.  PROGRAM  INFORMATION 
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TABLE  III.  SUBSONIC  SUBROUTINE 
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4.  SUPERSONIC  SUBROUTINES 


The  subroutine  that  directs  the  computation  for  the  supersonic  case 
is  shown  in  Table  IV.  It  differs  from  the  subsonic  subroutine  in  the 
computation  of  the  first  (potential)  term.  It  also  provides  values  for 
(dCN/da)  =  0  and  '(dC  /da)  =  0  when  they  have  not  been  supplied  as  data 
in  N  thea  main  program.  a 

5.  SECOND  TERM  INTEGRAL 

The  subroutine  that  computes  the  integral  portion  of  the  second 
(viscous)  terms  of  the  normal  force  and  moment  coefficients  is  shown  in 
Tabled.  It  initializes  incremented  variables,  sets  the  increment  for 
trapezoid  integration  for  the  nose  portion,  and  sets  the  initial  value  of 
x  as  one-half  an  increment  ahead  of  the  nose.  After  incrementing  x,  it 
computes  the  radius,  selects  the  drag  coefficient,  corrects  it  for  transient 
effect  and  finite  length,  then  computes  the  incremental  value.  When  the 
value  of  x  reaches  the  last  increment  on  the  nose  portion,  it  is  increased 
a  half- increment  to  equal  the  nose  length,  the  increment  size  is  changed, 
x  is  reset  one-half  increment  forward  of  the  beginning  of  the  afterbody, 
and  the  integration  is  continued.  For  the  afterbody  portion,  the  radius 
remains  constant  and  is  defined  outside  the  integration  loop. 

6.  CROSS-FLOW  DRAG  COEFFICIENT 

The  function  that  computes  the  drag  coefficient  for  the  integration 
subroutine  is  shown  in  Table  VI.  Using  the  cross-flow  Mach  number,  the 
function  selects  one  of  eight  equations.  Then,  if  the  cross-flow  Mach 
number  is  below  0.55,  the  free-stream  Reynolds  number  is  computed.  If  it 
is  in  the  transition  region,  the  drag  coefficient  is  adjusted  by  one  of 
two  functions. 

7.  TRANSIENT  EFFECTS 

The  function  that  corrects  the  cross-flow  drag  coefficient  for  the 
transient  effect  is  shown  in  Table  VII.  If  the  cross-flow- Mach  number  is 
greater  than  one,  no  correction  is  made.  Otherwise,  the  time  parameter  is 
calculated,  and  one  of  six  functions  is  selected  to  compute  the  correction 
factor. 

8.  RADIUS 

The  function  that  computes  the  radius  at  any  point  on  the  model  nose 
is  shown  in  Table  VIII.  The  proper  nose  type  is  determined  from  the  nose- 
type  code,  then  the  appropriate  function  is  selected.  For  the  power  series 
and  parabolic  noses,  a  simple  relation  is  used  to  generate  an  exponent,  and 
the  radius  is  computed.  The  ogives  are  represented  by  the  function  for 
tangent  ogives.  Radii  for  cones  are  computed  from  the  tangent  of  the  nose 
half-angle. 
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TABLE  IV.  SUPERSONIC  SUBROUTINE 
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TABLE  V.  SECOND  TERM  INTEGRAL  OF  NORMAL  FORCE  AND  MOMENT 
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TABLE  VI.  DRfiC-  COEFFICIENT  FOR  CYLINDER  AT  CROSSFLOW  MACH  AND  REYNOLDS  NUMBERS 
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********  ***********'*  ***************  *  *****************  ♦*  +  ♦*****♦**♦♦*** 


TABLE  VII.  CYLINDER  CROSS-FLOW  DRAG  COEFFICIENT  TRANSIENT  FACTOR 


ft*  — 


>0  0  0' 

a  o-  o 

o  u  *• 

H  ‘  K 

O  ^ 

O  O  s  u, 

in  in  <o  in  -K- 

o  -<  -i  -> 

H  •  •  O  • 


O  \ 
o  in/  — 
-1  a 

•  ►- 


>  s 
—>  o  ■— 
a  o  a 


->  a 
u  <  s 
o  — 
u  2  < 

w 

2  -J)  7. 
O  #  « 
—•  £< 

^  <  *  I 
U  It  X  ' 
Z  U  l| 

3  x  ai 

u.  <  H  ' 


•  i-  • 
-«  #  + 

sow 


II  II  II  II  -4  II  u 
II  u 
u .  u.  u. ,  u.  u.  u.  n 
U  U  U  U  U  u  u 

Q  O'/O  Q<Q  Q  Q 

yuuuuuy 

N  -*  <VI  M  <  (fl  «{>• 


25 


TABLE  VIII.  RADIUS  AT  ANY  GIVEN  X 


FUNCTION  RADIUS  (NtX.ALENN) 

IF  <N.£0. I )  GO  TO  2 
IF  (N.F-,0.2)  GO  TO  1 
P=ALOG(  .iS)/ALOG<  ALENN) 

RADIUS=X*»P  GO  TO  3 

1  0R=ALENN**2+.2S 

Y=SdRT<OR**2-(X-ALENN>**2> 
RAl)IUS-Y”QR+«S  GO  TO  3 

2  RAOIUS=< «5/ALENN>*X 

3  RETURN  END 

+  ****:***********♦**♦  ■♦♦♦*********’ 


9.  THETA 

The  theta  subroutine  that  computes  the  average  change  in  radius  with 
respect  to  axial  length,  is  shown  in  Table  IX.  First,  cones  are  considered, 
and  the  nose  ang,le  is  computed  directly.  Next,  incremented  variables  are 
initialized,  and  the  increment  is  computed.  Using  the  radius  function,  the 
radius  is  computed  at  each  increment  and  compared  to  that  of  the  previous 
increment .  These  values  are  summed  to  obtain  an  average  nose  angle.  Finally, 
the  nose  angle  is  proportioned  to  the  total  model  length  to  get  an  average 
angle  for  the  model. 

TABLE  IX.  THETA  FOR  THE  NOSE  MODEL  BEING  CONSIDERED 


SUBROUTINE  THETA  ( N  »  ALLNN  ,  ALENM  ,  -THEN.  THEM ) 

IF  (N.EQ.l)  GO  TO  2 

THEN  =0.  X=0*  D£LX=ALENN/10«  R0=0« 

1  X  =  X  +DELX 

RN=RA01US(N»X  »  ALENN ) 

RC=RN-RU 

THE  T=AT  AN (RC/DELX ) 

THEN  =THEN  +.I*THET  RO=RN 

IF  <  X«LT«  < ALENN-DEL X/2#  )  )  GO  TO  1  ..  10  3 

2  THEN  =ATAN< .5/ALENN) 

3  THEM  - 1 HEN  *  ALhNN/ALENM 

4  RETURN  END 

************************************************* 
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10.  DRAG  CORRECTION  FOR  FINITE  LENGTH 

The  function  that  corrects  the  two-dimensional  drag  coefficient  for 
finite  length  is  shown  in  Table  X.  The  length  from  the  nearest  end  is 
computed  first,  then,  the  cross-flow  Mach  number  and  the  square  root  of 
the  end  length  are  used  to  compute  a  correction.  If  the  point  being 
considered  is  oh  the  nose,  another  adjustment  is  made. 

TABLE  X.  DRAG  COEFFICIENT  CORRECTION  FOR  FINITE  LENGTH 


FUNCTION  CCCRR  (  X  ,  ALENM  ♦  THE  N  .  AW  ,  A  ,  At.CNN  ) 

FL=X  IF  tX.liT.ALf:NM/2.  )  KL=  ALCNM  -  X 
AMC=AM*S  IN.<  A  )  ARGf:XP=<  AMC*  1  ,  )*SORT{EL> 

OCORR=l  *-F.  AP  { —AK'GhXP  ) 

IF  I  X  «LT  •  ALCNN  )  UCCJRH  =  DCORR*CCJS(  THUN  ) 

RETURN  END 

*******  ♦  ♦♦♦♦****<^**  ***********  *♦***«■♦*♦*****# 
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11.  DCNDA 

The  DCNDA  function  that  computes  the  initial  slope  of  the  normal 
force  coefficient  when  plotted  against  the  angle  of  attack  is  shown  in 
Table  XL  After  the  afterbody  length  is  determined,  a  branch  to  the 
correct  Mach  number  is  taken.  For  each  Mach  number  branch,  a  function  is 
selected  to  correspond  to  nose1  type  and  afterbody  length.  Finally,  the 
value  is  converted  from  AC^  per  degree  to  AC^  per  radian. 

12.  DCMDA 

The  DCMDA  function  that  approximates  the  initial  slope  of  the  moment 
coefficient  when  plotted  against  angle  of  attack  is  shown  in  Table  XII. 

It  computes  the  length  of  the  moment  arm  from  a  function  which  describes 
the  movement  of  the  center  of  pressure  versus  Mach  number.  This  arm  is 
then  applied  to  the  DCNDA  supplied  from  the  argument  list  to  obtain  DCMDA. 
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TABLE  XI.  DCNDA 


rc  COMPUTE  OCNOA «  IF  NOT  SUPPLIED* 
THE  T A-SUD-V  FROM  5  TO  15  DEGREES 

FUNCTION  CNA-  < N , ALENN , ALENM, AM ) 

ALENA-=ALENM— ALfiNN 

IF  ( AM.L T .2.5  )  GO  TO  10 


IF  (AM.LT  .3.5)  GO  TO  6 
IF  ’(AM.LT.4.5)  GO  TO  2 
IF  ( ALENA . GT • A • )  GO  TO  1 
CNA-  .03?  +  . 00 J*ALCNA  - 

1  CNA=,043  GO  TO  14 

2  IF  ( ALENA .G  T .4.5)  GO  TC 
IF  (  N  .GT  .  I  )  GO  TU  3 

CNA=  .  0  3  J  ♦  .  004*ALt:  NA  ■ 

3  CNA  =  . 0  J  8  ♦  .00  J* ALENA  - 

4  IF  ( ALENA. GT .H.O )  GO  TO 
CNA=  .052  +  .001 *ALENA  - 

5  CNA=.047  GO  TO  14 

6  IF  I  ALENA. GT  .3.5)  GO  TO 
IF  ( N  .GT . 1 )  GO  TO  7 
CN4=  .033  +  •  0 0.6 GALENA  - 

7  CNA=  .043  ♦  • 0044ALLNA  - 
0  IF  ( ALENA.  GT .8.0)  GO  TO 

CNA=  .056  «■  .001 *ALLNA  - 
9  CNA=.049  GO  TU  14 
10  IF  (  ALENA  .vjT  .2.5)  GO  TO 
IF  (N.GT.l)  GO  TO  ll 
CNA-  .033  4-  •  0  l  04 ALENA  - 
■  11  CNA=  .043  ♦  .0054ALENA  - 

12  IF  ( ALENA .GT . T  .0  )  GO  TO 
CNA=  .063  +  • G0?*ALENA  - 

13  CN A= .045 

14  CNA=CNA*tfc0./3. 14159 

*********** **********4 

TABLE  XII. 

TO  COMPUTE  DCMDA 


• 02* ( Am— 5* )  GO  TO  14 

4 

•  04  * ( AM— 4 • )  CC  TO  14 

.  Ob  *  (  AM— 4  •  )  GO  TO  14 

5 

. 044 { AM— 4 . )  GO  TU  14 

8 

•  03*( AM-3. >  GO  TO  14 

.  04  * ( AM— 3  .  )  GC  TO  14 

9 

»  03* ( AM— 3 .  )  GO  TO  14 

12 

.03* ( AM— 2 . )  GO  TO  14 

.  04  * ( A M—  2  •  )  GO  TO  14 

13 

»  08* ( Am- 2  .  )  GO  TO  14 

RETURN  END 

*^**♦******♦♦♦****♦**4:*  ****** 

DCMDA 

.  IF  NOT  SUPPLIED 


FUNCTION  CMA  ( DCNDA , ALENM . AM ) 

ARM= ( . 75- . 0634AM )* ALENM 
CMA-DCNDA4ARM  RFTURN  END 

**  *!******  ******************  ******; 
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13.  OUTPUT 

The  output  format  is  shown  in  Tables  XIII  and  XIV.  The  values 
printed  at .zero  angle  of  attack  are  to  remind  the  reader  that  the  procedures 
used  are  ail  based  on  simple  bodies  of  revolution  which  create  only  drag 
(no  normal  force  or  moment)  when  the  angle  of  attack  is  zero.  The  drag 
coefficient  printed  is  that  of  the  last  increment  of  integration  on  the 
model  afterbody,  without  corrections  for  transient  effect  and  finite  length 
effect., 


TABLE  XIII.  OUTPUT  (LEFT  HALF) 


HUN  OF  NOSE  CODt£  1  MODEL.  WITH  A  2.83  CALIBER  NOSE. 

AND  AN  OVERALL  LENGTH  OF  O.M.i  CALIBERS.  MACH  NUMBER  IS  6.04 
AND  Rll  ( PER  CAL)  OVER  MACH  NO.  IS  22700. 


ALPHA 

CSUBN 

CSUBM 

FSTTMN 

DCNOA  IS 

2 . 46  »  AND  DCMDA 

IS  8.43 

0.00 

0.  00 

0.00 

5.0  0 

^.24 

Li.  8.1 

0.21 

10.00 

0.85 

3.11 

0.42 

15.00 

1.43 

5.22 

0.62 

2  0.00 

2.  14 

7.84 

0.79 

25.00 

2.97 

10.93 

0.94 

.10.00 

3.92 

14.44 

1.07 

35.00 

4.92 

18.17 

1.16 

40.00 

5.95 

22.00 

1.21 

45.00 

6.97 

25.8  2 

1  .23 

50.00 

7.95 

29.51 

1  .2! 

55.00 

CJ.87 

32.96 

1.16 

60.00 

9.69 

36.07 

1.07 

65.00 

10.39 

30.73 

0.94 

70.00 

10.95 

40.88 

0.79 

75.00 

11.35 

42.44 

0.62 

80.00 

11.58 

43.37 

0.42 

85.00 

11.63 

43.63 

0.21 

90.00 

11.50 

43.23 

0.00 
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TABLE  XIV.  OUTPUT  (RIGHT  HALF) 


SECTMN 

FSTTKM 

0.02 

0.73 

0.43 

1.44 

0.81 

2.11 

1.35 

2.71 

2.03 

3.23 

2.85 

3.65 

3.76 

3.96 

4.74 

4.15 

5.74 

4.21 

6.74 

V.15 

7.71 

3.96 

8.62 

3.65 

9.45 

3.23 

10.16 

2.71 

10.73 

2.11 

11.16 

1.44 

1  1  .42 

0.73 

1  1  .50 

0.00 

SECTMM 

CO 

0.08 

1.59 

1.67 

1.59 

3.11 

1.32 

5.13 

1.24 

7.70 

1.22 

10.79 

1.22 

14.21 

1.22 

17.85 

1.22 

21.61 

1.22 

25.36 

1.22 

29.00 

1.22 

32.42 

1.22 

35.50 

1.22 

38.17 

1.22 

40.33 

1.22 

41.92 

1.22 

42.90 

1.22 

43.23 

1.22 

SECTION  V, 
RESULTS 


The  computer  program  was  executed  more  than  40  times,  with  over  25 
runs  producing  useable  data.  Many  paired  runs  .were  made  to  observe  the 
effect  of  changing  a  single  parameter  or  constant.  Near  the  final  stages 
of  the  program  development,  a  separate  regression  program  was  iL  4  to  fit 
one  set  of  experimental  data.  A  function  of  the  form  (constants)’  sin  a 
cos  a  +  (constant)  sin3a/cos  a  was  used  for  this  regression  program, 
since  it  represents  the  alpha -dependent  quantities  of  the  theory.  The 
high  degree,  of  conformity  indicated  the  method  was  of  the  proper  form. 

Three  sources  were  used  for  comparisons  with  the  output  of  the  program. 
Almost  no  data  for  the  desired  angle-of-attack  range  are  available. 

Additional  data  are  required  for  angle  of  attack  up  to  90  degrees,  parti¬ 
cularly  at  transonic  speeds. 

The  mpment  coefficient  will  not  be  discussed  in  detail,  since  the 
normal  force  is  calculated  for  incremental  areas,  and  if  it  is  accurate 
and  no  couples  exist,  the  moment  must  also  be  accurate.  This  can  be 
observed  from  comparisons  with  the  moment  data  in  References  15  and  16, 
where  deviations  in  normal  force  coefficient  are  reflected  almost  directly 
in  the  moment  coefficient. 

1.  SUBSONIC 

The  subsonic  runs  were  compared  with  Hauer  and  Kelly's  data  in  Reference 
15  which  reports  primarily  on  the  Magnus  effect  of  spinning  projectile 
models,  but  also  presents  data  for  a  zero  spin  rate.  As  shown  in  Figure  11, 
the  computed  values  of  the  coefficients  generally  take  the  same  shape  as 
the  experimental  data  up  to  about  25  degrees.  From  there  to  about  55 
degrees  ,  the  data  show  a  marked  divergence.  Hauer  and  Kelly  describe  a 
region  of  flow  unsteadiness  at  55  degrees  which  is  attributed  to  a  rapid 
transition  between  an  axial  and  a  transverse  flow  pattern.  Above  this 
region  of  unsteady  flow,  the  comparison  shows  another  divergence,  but  at 
a  different  degree  of  error. 

Much  of  the  error  in  the  comparison  seems  to  come  from  the  values  of 
the  cross-flow  drag  coefficient  used  since  the  computations  gave  excellent 
comparisons  to  55  degrees  and  then  severe  overpredi  colons  to 90  degrees,  if 
a  laminar  value  were  used.  The  difficulty  lies  in  the  uncertainty  involved 
in  predicting  the  transition  from  laminar  to  turbulent  flow,  and  the  rate 
of  recovery  of  the  drag  coefficient  above  the  transition  point.  Reference 
11  presents  several  different  cases  which  could  apply,  but  the  selection 
of  the  proper  one  would'  probably  be  based  on  prior  knowledge  of  the 
experimental  data.  Transition  for  a  cylinder  in  cross-flow  usually  occurs 
at  a  diameter-based  Reynolds  number.^  500,000,  but  this  may  vary  by  as 
much  as  200,000.  If  the  test  case  lies  within  this  region  of  Reynolds 
numbers,  the  actual  transition  Reynolds  number  would  have  to  be  known  to 
permit  selection  of  the  proper  cross-flow  drag  coefficient. 
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Figure  11.  Comparisons  to  Experimental  Data-Subsonic  (Reference  15) 
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Figure  11.  (Concluded) 
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2.  SUPERSONIC 

The  supersonic  runs  (Figure  12)  were  made  with  an  input  value  for 
DCNDA.  The  slope  used  for  input  either  was  taken  from  the  curve  in 
Reference  14,  or  was  computed  from1  the  experimental  data  between  zero  and 
five  degrees  angle  of  attack.  This  computation  can  be  used  because  the 
various  runs  show  that  potential  effects  account  for  95  percent  or  more 
of  the  total  normal  force  coefficient  in  this  region. 

The  runs  show  excellent  agreement  with, the  data  presented  by  Penland 
in  Reference  13.  The  agreement  in- the  hypersonic  region  is  excellent, 
since  the  flow  approaches  the 'conditions  that  satisfy  the  Newtonian  theory. 
The  agreement  is  slightly  better  than  that  produced  by  the  Modified* 
Newtonian  theory.  The  data  used  were  in  the*  form  of  lift  and  drag  coeffi¬ 
cients  referenced  to  planform  area,  but  a  short  program  on  desk-topv 
computer  converted  the  data  to  normal  force  coefficient  referenced  to  base 
area.  The  error  introduced  by  these  computations  is  not  significant  for 
the  purpose  of  comparison. 

The  comparisons  at  Mach  2.37,  2.98,  and  3.90  require  additional  con¬ 
sideration.  In  an  effort  to  determine  n»  the  90* degree  angle-of-attack 
runs  with  cylindrical  noses  were  used  to  ascertain  the  effects  of *Mach 
number.  When  the  experimental  data  were  converted  useable  form,  the 
drag  coefficient  exceeded  the  value  for  an  infinite  cylinder.  Such  a 
condition  is  contrary  to  all  previous  knowledge  of  this  ttype  of  flow,  and 
an  error  is  indicated.  Since  the  test  runs  Were  made  with  a  wall-mounted 
sting,  the  corrections  for  the  mounting  condition  may  contain  some  error. 
However,  the  data  are  useful  for  plotting  the  shape  of  the  curve  and  were 
adequate  for  the  original  rcoort.  The  wall  mounting  limited  the  data 
start  to  30  degrees  and  prevented  a  check  on  the  test  apparatus,  at  zero 
angle  of  attack.  Therefore,  the  source  of  the  apparent  error  cannot  be 
evaluated. 

3.  INDICATED  IMPROVEMENTS 

As  the  program  developed,  several  areas  were  indicated  where' an 
increase  in  complexity  might  improve  the  results;  however,  the  return- fromi 
the  additional  effort  would  be  low.  Also,  the  lack  of  experimental  data 
for  correlation  with  the  theoretical  results  precluded:  any  real  value  in 
an  attempt  to  improve  the  accuracy* of  the  method.  If  valid  data  were 
available  and  if  greater  accuracy  were  required,  some;  improvement  in  the 
results  would  be  possible. 

The  best  candidate  area  for  improving  accuracy  would.be  the  drag 
coefficient  in  the  region  of  transition  Reynolds  numbers.  A  method  of 
reading  into  the  computer  the  transition  Reynolds  number,  based  on  surface 
roughness  and  other  pertinent  factors,  could  be  used  to  adjust  the 
transition  induced  dip  in  the  C .  versus  Reynolds  number  curve.  Certainly, 
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★  Experimental  □  Computed 

♦  Modified  Newtonian  Theory 
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Figure  12.  (Continued)  (Reference  17) 
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additional  experimental  data  in  the  subsonic  range  would  be  required  for 
a  significant  improvement  in  this  function.  More  attention  to  the  change 
from  axial  to  transverse  flow  dominance  could  lead  to  a  factor  to  define 
its  manner  and  causes.  Since  both  the  axial  and  the  transverse  flows 
are  involved  in  transition,  and  since  transition  is,  at  best,  not  a  very 
predictable  phenomenon,  a  reliable  prediction  may  not  be  possible. 

There  are  other  improvement  areas  which  would  give  lesser  benefits  in 
the  improvement  of  accuracy,  and  they  are  mentioned  in  the  order  in  which 
they  appear  in  the  program. 

The  volume  in  the  subsonic  subroutine  assumes  an  inscribed  cone  for 
noses  which  are  not  true  cones.  The  radius  function  could  be  expanded, 
or  an  additional  function  could  be  added,  to  provide  a  value  of  volume 
that  would  be  more  accurate  for  non-conical  noses.  In  particular,  this 
modification  would  yield  better  results  for  the  blunter  class  of  power 
series  and  parabolic  hoses. 

Since  the  DCNDA  function  is  approximate,  it  could  easily  be  improved 
with  additional  experimental  data.  The  insertion  of  experimental  data 
on  the  input  data  cards  would  be  the  easiest  method  of  correcting  this 
source  of  error. 

The  theta  subroutine  is  an  averaging  process,  and  it  is  necessary  that 
theta  be  averaged  for  the  apparent  mass  factor.  However,  when  using  the 
nose  angle  to  correct  n,  the  local  value  of  theta  would  be  more  accurate 
for  noses*  which  are  not  true  cones. 

All  the  functions  that  are  approximated  with  line  segments  could  be 
fitted  with  higher  order  curves  that  would  be  more  accurate. 

Any  of  the  factors  mentioned  might  produce  some  refinement,  in  the 
accuracy  of  the  predicted  aerodynamic  characteristics;  however,  the  added 
complexity  is  not  justified  at  this  time,  since  experimental  data  for 
actual  verification  are  not  available.  In  addition,  some  of  the  recommend¬ 
ed  techniques  require  experimental  data  for  their  implementation.  Further, 
it  is  a  primary  goal  of  this  study  to  retain  a  minimum  of  complexity. 
Additional  study  should  be  performed  to  improve  the  technique,  but  without 
significantly  increasing  the  level  of  complexity. 
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SECTION  VI 


FUTURE  DEVELOPMENT 


The  most  significant  improvement  that  might  be  made  to  the  program 
would  be  to  include  more  complex  models.  The  first  additions  might  be 
nose  blunting  and  boattailing  of  the  afterbody,  since  each  of  these  could 
be  quite  readily  incorporated  in  the  present  program  with  only  minor 
modifications.  The  addition  of  fins,  wings,  control  surfaces,  model 
irregularities,  and  changes  in  shape  to  include  other  than  bodies  of 
revolution  are  some  of  the  extensions  that  would  be  desirable.  These 
additions  require  extensive  effort  for  theoretical  development,  as  well 
as  large  amount  of  experimental  data;  however,  the  results  would  be  most 
valuable  to  future  designers  and  analysts. 

The  requirements  of  the  Air  Force,  the  intended  use  of  the  program, 
and  the  success  of  early  advances  in  complexity  would  determine  how  much 
extension  should  be  undertaken.  At  this  time-*  the  program  development 
confirms  the  basdc  theories  involved  and  verifies  the  value  of  the  com-  ■ 
puter  as  a  basic  research  tool  in  this  area. 
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APPENDIX 


The  following  block  diagrams  are  intended  to  show  the  various  'relation¬ 
ships  of  the  different  parts  of  the  program  to  each  other  and  to  thb  Main 
Program.  Only  those,  subprograms  which  utilize  another  subprogram  are 
included.  These  diagrams  are  generalized  by  the  various  parts  of  the  sub¬ 
program  being  considered  and  should  not  be  construed;  to  be  flow  charts. 

The  equations  necessary  to  write  the  flow  charts  have  already  been  presented. 


Figure  1-1.  Main  Program 
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Figure  1-2.  Subson  Subprogram 
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Figure  1-3.  Supson  Subprogram 
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Figure  1-4.  Secti  Subprogram 
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